Itinerant ferromagnetism in the multiorbital Hubbard model 
a dynamical mean-field study 
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In order to resolve the long-standing issue of how the itinerant ferromagnetism is affected by the 
lattice structure and Hund's coupling, we have compared various three-dimensional lattice struc- 
tures in the single- and multiorbital Hubbard models with the dynamical mean-field theory with 
an improved quantum Monte Carlo algorithm that preserves the spin-SU(2) symmetry. The result 
indicates that both the lattice structure and the d-orbital degeneracy are essential for the ferromag- 
netism in the parameter region representing a transition metal. Specifically, (a) Hund's coupling, 
despite the common belief, is important, which is here identified to come from particle-hole scat- 
terings, and (b) the ferromagnetism is a correlation effect (outside the Stoner picture) as indicated 
from the band-filling dependence. 

PACS numbers: 71.10.Fd; 75.20.En; 75.40.Mg 



Itinerant ferromagnetism in transition metals and 
their compounds has been one of the central issues in 
condensed-matter physics. Despite a long history of in- 
vestigations, dating back to e.g., the Stoner 1 or Slater- 
Pauling 2 theories, a full consensus on the mechanism has 
not been achieved yet. While the single-orbital Hubbard 
model was introduced as the simplest model to capture 
the ferromagnetism from electron-electron interactions, 
it has become increasingly clear that the model on usual 
lattices does not show ferromagnetism for realistic val- 
ues of the Coulomb repulsion U. Intensive studies have 
then ensued to incorporate the ingredients other than the 
Hubbard U . Two factors are now considered to be sig- 
nificant: lattice structure and the degeneracy of d orbits. 

Early theories due to Kanamori 3 and Gutzwiller al- 
ready suggested that lattice structure (which dictates the 
shape of the density of states (DOS)), is crucial for the 
ferromagnetism. Specifically, they discussed the itiner- 
ant ferromagnetism for a face-centered cubic (fee) lattice 
as in Ni where the DOS is peaked at the band edge, 
and argued that the peak stabilizes the ferromagnetism. 
The problem is revisited in recent studies: When the 
orbital degeneracy is ignored, Arita et al. showed with 
the fluctuation exchange (FLEX) and the two-particle 
self-consistent approximations that a ferromagnetic ten- 
dency is strongest for fee. The ferromagnetism on fee lat- 
tices has also been found by Ulmke 6 with the dynamical 
mean- field theory (DMFT). 7 The importance of lattice 
structure has also been suggested by some exact results 
for flat-band systems. 8 

On the other hand, the importance of the cZ-orbital 
degeneracy and the associated Hund exchange coupling 
has long been stressed. 9 While Hund's coupling is an 
intra-atomic interaction favoring aligned electron spins, 
the interaction may cause a long-range ferromagnetic or- 
der through electron transfers. This was followed by 
intensive studies, but the complexity of a multiorbital 
model has limited the existing studies to some restric- 
tions: The two-orbital Hubbard model in one dimension 



(ID) has been most intensively studied, where a ferro- 
magnetic ground state with an antiferro-orbital order is 
expected for the quarter filling (n = 1; one electron per 
site), for strong repulsion U and Hund's coupling J. In- 
deed, quantum Monte Carlo (QMC) 10 , exact diagonal- 
ization (ED), 11-13 and density-matrix renormalization- 
group (DMRG) 14 studies for finite-size chains have con- 
firmed the ferromagnetic ground state. However, these 
quarter-filled systems are insulating. The purpose of the 
present letter is to look into metallic ferromagnetism in 
multiorbital systems, with our eyes set on transition- 
metal ferromagnets such as Ni. 

Several studies on metallic ferromagnetism in multi- 
orbital systems exist. 11,14-18 Numerically exact results 
for the double-orbital Hubbard model in ID have been 
obtained with ED 11 - 13 or DMRG 14 , while those on 
the infinite-dimensional hypercubic 16 or Bethe 17 lattices 
have been studied with the DMFT. These studies have 
found itinerant ferromagnetism away from the quarter 
filling only for very large Hund's couplings J > W/2 (W: 
bandwidth), and no ferromagnetism has been found in 
a realistic range of J for transition metals. However, 
it is in our view still an open question whether Hund's 
coupling is essential in real ferromagnets, since these cal- 
culations do not take account of lattice structures of real 
materials. 18 

This has motivated us to investigate here the effect 
of lattice structures and of multiorbital correlations with 
the DMFT, comparing a simple-cubic and fee lattices in 
the single- and multiorbital cases. We shall conclude that 
both the lattice structure and Hund's coupling are crucial 
for itinerant ferromagnetism in realistic (~ transition- 
metal) parameter regions. Physically, we observe that 
(i) particle-hole scatterings, neglected in Kanamori's T- 
matrix theory, are essential in fact in the presence of 
Hund's coupling, and (ii) the present result for the band- 
filling dependence of the magnetism is totally out of 
Stoner's picture. 

So we start from the multiorbital Hubbard model on 
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FIG. 1: (Color online) The density of states for noninteracting 
electrons on the fee lattice. Red figures represent the band 
filling (per spin), while black ones the corresponding total 
band filling for a two-orbital system. 



three-dimensional lattices. For the fee lattice the disper- 
sion is 

3 

e(k) =4t^cos(fci)cos(fcj) + 2£'^cos(2fci), (1) 

where t(t') is the nearest-neighbor (second-neighbor) 
hopping, and a cubic Brillouin zone (— n < ki < it) for 
two equivalent, interpenetrating fee lattices is adopted. 
We set t = 4t' = -^7= to fix the effective bandwidth at 

3vll 

W = 4(J e 2 D(e)de) 1 / 2 = 4. For simplicity we have as- 
sumed isotropic transfer integrals and ignored hybridiza- 
tions between orbitals. The DOS for noninteracting elec- 
trons on fee in Fig. 1 has a peak at the lower band edge, 
where Ep resides around the peak in Ni in the hole pic- 
ture. For comparison we also consider the simple-cubic 
lattice, for which the dispersion is e(fc) = 2t^ j=1 cos(/ci) 
with t = -4 and W = 4. 



The interaction Hamiltonian is 

-ffint = Hjj + Hj, 



(2) 



i,m<m f ,a 

Hj — J ( c imT C im' l c imlCim't + c l m f c Lil C ""'l C iTO'T )' 

where c\ mcr creates an electron with spin a in orbital m at 
lattice site i, rii ma = cl ma Ci ma , U(U') denotes the intra- 
(inter-)orbital Coulomb interaction, and J the Hund- 
exchange and pair-hopping interactions. Here we have 
decomposed H m t into Hu (the density-density interac- 
tions) and Hj (not expressible as such). 

We investigate the above model with the DMFT, which 
neglects the fe-dependence in the self-energy and in the 
vertex but incorporates their temporal dependence and 
the one-electron dispersions. Since the DMFT becomes 
exact in the limit of large coordination numbers, the ap- 
proximation is expected to be fair for fee, where a site has 



twelve nearest neighbors and six second neighbors. The 
DMFT impurity problem is solved exploiting the QMC 
method developed in our previous paper. 19 The QMC 
algorithm combines the Trotter decomposition for Hjj 
and a series expansion for Hj to decouple the two-body 
interactions with the auxiliary-field transformations. 20 
While the algorithm uses a series expansion, it is vir- 
tually nonperturbative, since all the nonvanishing orders 
are incorporated numerically. An important virtue of 
the present QMC is that it preserves the spin (SU(2)) 
and orbital (U = U' + 2 J) rotational symmetries in Hi nt , 
which arc difficult to preserve in the conventional QMC 
method. 17 ' 21 The present method also enables us to ad- 
dress temperatures close to the Curie temperature Tc, 
to which the ED or the DMRG cannot access. 
We calculate the spin susceptibility, 



X(0,0) = ]Txw(O,0), 
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in the paramagnetic phase (T > Tc) in the two-orbital 
Hubbard model, through the Bethe-Salpeter equation, 



(3) 



where x IS the susceptibility matrix with respect to 
to, to', iuj 1 , xo the irreducible lattice Green function 
and f the vertex function. The DMFT approximates 
the electron self-energy S(fe, iu>) and the vertex function 
T(k,iuj;k' ,iuj') to be local, T,(iuj) and T(iuj,iuj') respec- 
tively, which are obtained from a QMC calculation for 
the self-consistently-determined impurity model. 

First, we discuss the results for the simple-cubic lat- 
tice (not shown). We used (U,U',J) = (8,6,1), which 
are reasonable (or somewhat overestimated) values for 
transition metals. x(®> 0) does not exhibit divergent be- 
haviors at any filling even when we extrapolate the results 
to low temperatures, which indicates that the itinerant 
ferromagnetism cannot be solely attributed to the mul- 
tiorbital interactions. This result is similar to the result 
for an ooD hypercubic lattice, 16 where the ferromagnetic 
ground state appears only for very large U > 3W and 
J > W/2. 

Now we turn to the fee lattice. Figure 2 shows how 
the temperature dependence of x(0, 0) _1 for n = 1.5, 0.4 
changes as we successively introduce U' [(U, U', J) = 
(4,0,0): single-orbital case — » (4,4,0)], and then J 
[(4,4,0) — > (4, U', J)], preserving the rotational symme- 
try (U = U' + 2 J) for multiorbital cases (W ^0). We 
can see that x is noticeably suppressed and the transi- 
tion disappears when U' alone is switched on. When J 
is added, however, x is significantly enhanced, and Tc 
is pulled back to finite values comparable to that in the 
single-orbital case. The way in which Tc is pulled back 
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FIG. 2: (Color online) Inverse of the spin susceptibility x(0, 0) 
against temperature for the two-orbital Hubbard model on 
the fee lattice for U — 4 for band filling n = 1.5(a) and 
n = 0.4(b). We switch on U (blue) and then J (orange 
and red) preserving the rotational symmetry (U = U' + 2J). 
The solid lines are guide to the eye, and the dashed lines 
extrapolations. 



strongly depends on J: At n = 1.5, Tc — 0.015 for 
J = 0.25 while T c ~ 0.05 for J = 0.5. The result clearly 
shows a crucial role of the interorbital interactions (U' 
and J) in the ferromagnetism. 

In Ni the triply-degenerate t^g bands have 0.6 holes per 
site, so that n = 0.4 for the present two-orbital model 
roughly corresponds to the band filling of Ni in the hole 
picture. Thus the present result, Fig. 2(b), implies a sub- 
stantial role of Hund's coupling in the metallic ferromag- 
netism in the parameter region that contains Ni's. This 
is to be contrasted to the Kanamori theory, to which we 
shall come back. 

Figure 3 displays the band-filling dependence of x(0, 0) 
at T = 0.1 for various values of U = 2 — 5, where we have 
fixed J = 0.5 and the relation U 1 = U — 2 J. The sus- 
ceptibility is seen to take the largest value at around the 
filling rip ~ 1.2-1.5, where n p shifts markedly to higher n 
as U is increased. So the filling dependence of x is very 
dissimilar to the initial DOS (Fig. 1), which has a peak at 
around n ~ 0.4. This is to be contrasted with one-band 
results such as Stoner's, T-matrix 3 , or FLEX theories. 5 ' 22 
The result is even distinct from multiorbital results such 
as the Gutzwiller approximation, 18 where a higher DOS 
at Ep favors ferromagnetism. Hence the present result 
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FIG. 3: (Color online) Spin susceptibility against the band 
filling n for the two-orbital Hubbard model on the fee lattice 
at T = 0.1 for various values of U, with J = 0.5 and the 
relation U' = U — 2J fixed. 
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FIG. 4: (Color online) Wave number dependence of the sus- 
ceptibility on the fee lattice for various values of n — 1.25 — 2.0 
for (7 = 4, U' = 3, J = 0.5, T = 0.1. 



indicates that the correlation effect is indeed involved in 
the ferromagnetism. 

Here we identify a component that contributes to the 
correlation effect. In the Kanamori theory with the T- 
matrix approximation, which takes only particle-particle 
(p-p) scatterings, the direct interaction is reduced to 
17/(1 + U X PP ) while Hund's term to [J/(l + U X PP ) 2 }, 
so the latter is concluded to be less relevant. However, 
while the T-matrix approximation is justified only for low 
electron densities, other types of scatterings such as the 
particle-hole (p-h) channel 5,23 can play a role for general 
band fillings. In particular, the p-h scatterings enhance 
U into U/(l- Ux ph ) (x ph; the p-h propagator), and J 
is even more enhanced into [J /(I - Ux ph ) 2 }- 

Since the 

present DMFT calculation effectively includes these p-h 
scatterings, this should be one component that substan- 
tially contributes to the ferromagnetism. 

The interaction JJi nt should become more effective, es- 
pecially in the p-h channel, as the band filling approaches 
the half filling (n — 2). This will enhance x, as actually 
seen in Fig. 3 for n < n p . n p shifts to higher densities 
with U. This should be because the enhancement of x 
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FIG. 5: (Color online) Local spin moment (black) versus J 
in the two-orbital Hubbard model on the fee lattice for (a) 
n = 1.5 and (b) n = 0.75 with T = 0.1, U = 4 and the 
relation U' = U — 2 J fixed. The corresponding susceptibility 
x(0, 0) is plotted in blue. 

with U becomes stronger toward the half filling, while the 
original DOS is conversely peaked around low n. On the 
other hand, x starts to decrease for n > n p , which is ex- 
pected to come from antiferromagnetic correlations. To 
confirm the antiferromagnetic correlations, we have cal- 
culated the wave number dependence of 0), where 
the q dependence incorporated through xo in Eq. (3). 
We can see that \ m Fig- 4 has a peak at around T 
point for n = 1.25 and 1.5, while the peak disappears for 
n = 1.75 in favor of a diverging peak at around X and 
W points for n — 2, which indicates strong antiferromag- 
netic correlations as the half filling is approached. 

We make another argument from the the strong- 
coupling limit, where the Hamiltonian is given only 
by the intra-atomic interaction .Hint- The two-electron 
eigenstates are then classified by the spin and orbital 
symmetries as 



Notation Spin Orbital Expression Energy 

IS 1 * 1 singlet sym. 7^( c i-f c ij. ^ ^T ^!.) U±J 

IS singlet sym. tJtj ( c it c 2j. + c 2T c u) U' + J 

3°A triplet antisym. -^(cj^c^ — c^ii) U' — J 

3°A triplet antisym. c ict c 2<t U' — J. 

While all the spin-triplet states have the same energy 
V — J, the orbital-symmetric states split into three ener- 
gies, U+ J, U' + J and U ~ J. Then the ground state for 
(U, U', J) = (4, 0, 0) is a superposition of the four states, 
IS , and 3°' Ta A, while that for (4,4,0) becomes a su- 
perposition of the six states, 1S 0,± and 3°'^A. Namely, 
the latter has three times greater number of spin singlets 
than the former. This intuitively accounts for the reduc- 
tion of x(0, 0) with U' in Fig. 2. 

In the strong-coupling limit, J further brings down the 
energy of the spin triplets. This may be one reason for 
the increase of x(0, 0) with J in Fig. 2: The local spin 
moment, in Fig. 5 actually increases with J for n = 
1.5. For n — 0.75, however, the local moment remains 
almost constant while %(0, 0) noticeably increases with J. 
This implies that Hund's coupling, despite being a local 
interaction, aligns spins nonlocally via the one-electron 
hoppings. 

Future problems include the effect of anisotropic trans- 
fers or hybridization of d orbits, and an examination of 
the chemical trend in the itinerant ferromagnetism. 
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